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SENSITIVITY  ASSESSMENTS  OF  SOME  DISTURBED 
E AND  F REGION  REACTION  RATES 

I Introduction 


Tlie  disturbed  E and  F regions  of  the  ionosphere  play  important  roles 
on  systems,  detectors  and  communications.  To  assess  the  impact  of  such 
an  environment,  a detailed  understanding  of  the  physics  and  the  chemistry 
of  these  regions  is  essential.  .An  important  aspect  of  this  problem  is  the 
ionization  and  the  deionization  of  the  regions.  These  depend  critically 
on  a large  number  of  reactions  that  occur  in  the  disturbed  region  which 
in  turn  control  and  determine  the  extent,  duration  and  decay  of  the 
disturbance.  The  ionization  and  the  deionization  is  generally  described 
by  a time  dependent  code  which  monitors  the  density  distributions  of  the 
disturbed  ionospheric  constituents.  The  complexity  of  such  codes  depend 
heavily  on  their  utility  in  conjunction  with  large  scale  hydrodynamic 
codes,  the  desired  output  and  the  speed  of  calculation.  Two  NRL  Codes, 

The  Master  Code  and  the  Simple  Code  have  been  developed^  to  assess  the 
ionization,  deionization  and  radiation  from  the  E and  F regions. 

The  more  detailed  code,  the  .AIRE  Master  Code,  remains  dynamic  in  the 

sense  of  utilizing  up-to-date  reaction  rates.  It  can  be  utilized  to  assess 

the  impact  of  various  reaction  rates  on  the  ionization  and  the  deionization 

of  the  E and  the  F regions.  In  this  report,  we  utilize  the  current  .NRL 
2 

Master  Code  for  the  assessment  of  several  reaction  rates  on  the  deionization 
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N't,  0^,  02(a‘^7r),  NO"",  O'",  o''(^D),  0'"(^P),  N"",  N'"(^D)  and/c^S).  In 
addition,  the  Master  Code  calculates  three  temperatures;  the  electron 
temperature,  the  N,  vibrational  temperature  and  the  heavy  particle 
temperature.  The  time  histories  of  these  species,  the  electron  density 
and  the  relevant  temperatures  provide  the  ingredients  for  the  calculations 
of  the  radiation  emitted  from  the  disturbed  atmosphere. 

Ill  Results 


Six  sets  of  initial  conditions  were  chosen  and  the  current  NTIL 

7 

Master  Code"  was  utilized  to  carry  out  the  deionization  calculations. 
Figure  1 shows  the  time  histories  of  the  electron  densities  for  these  six 
cases.  Selecting  a special  case,  e.g.  Case  I,  the  time  histories  of  the 
atomic  ions  (Fig.  2),  the  molecular  ions  (Fig.  3),  the  atomic  metastables 
(Fig.  4),  and  the  densities  of  O^(a^A),  NO  and  the  three  temperatures 
(Fig.  4)  are  shown  in  Figures  2-5. 

In  keeping  a general  code  dynamic,  one  must  incorporate  most  current 
and  reliable  reaction  rates.  This  procedure  has  been  applied  continuously 
to  NRL  Master  Code  from  its  inception^.  Many  of  the  reaction  rates 
relevant  to  the  E and  F regions  of  the  ionosphere  and  are  utilized  in  NRL 
Master  Code,  have  been  well  established  under  the  reaction  rate  support 

4 

program  of  the  Defense  Nuclear  .Agency.  Some  of  the  reaction  rates, 
however,  require  further  experiments  so  that  their  values  will  be  well 
established  over  temperature  ranges  of  interest. 

In  this  report  we  assess  the  impact  of  certain  reaction  rates  on  the 
electron  density.  For  e.xample,  consider  the  dissociative  recombination 
no""  + e ^ .V  + 0 (1) 
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where  the  rate  coefficient  utilized  previously  in  the  .\RL  Master  Code  had 

-12  5 

a temperature  dependence  of  . This  value  was  obtained  by  Biondi 

fromseveral  measurements^’ However,  recently  two  new  measurements^ 

of  reaction  (1)  have  been  made  in  two  different  experiments.  These  two 

experiments  give  two  entirely  different  temperature  dependences  for  the 

9 -0  37 

rate.  Huang,  et  al  , give  a temperature  dependence  of  T^  ‘ while  Walls 

“0  83 

and  Dunn  obtain  a temperature  dependence  of  (T^)  ' . These  two  rate 

coefficients  are  shown  in  Figure  6 along  with  the  old^  rate  coefficient. 

It  is  seen  from  Figure  6 that  the  rate  coefficient  of  Huang,  et  al^  is 
large  compared  with  that  measured  by  Walls  and  Dunn^^,  even  though  the 
last  measurement  is  undoubtedly  for  .NO^  being  in  its  ground  vibrational 
state  {y=0).  The  rate  coefficient  with  its  dependence  on  both  the  electron 
and  the  vibrational  tenperature  should  follow 

(2) 

which  makes  the  two  rates  somewhat  of  a puzzle.  However,  recent  theoretical 
calculations  by  Lee^^,  attributes  the  discrepancy  as  being  due  to  resonances 
in  the  dissociative  recombination  cross  section. 

Using  the  three  rate  coefficients  for  the  dissociative  recombination 
of  NO  shown  in  Figure  6,  sensitivity  calculations  were  made  using  the  .NRL 
Master  Code  for  three  different  initial  conditions  indicated  as  Cases  I, 

IV  and  V.  The  results  for  the  electron  and  NO^  densities  are  shoAvn  in 
Figures  7-12.  It  can  be  stated  clearly  that  the  impact  of  this  rate 
coefficient  is  very  important  on  the  deionization  of  the  disturbed 
atmosphere,  and  as  such  on  detection  systems  and  communications. 

In  a highly  disturbed  E and  F region,  the  concentrations  of  excited 


Coefficient  ~ T (1 
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metastables  are  appreciable, and  a large  number  of  near  resonance  charge 

exchange  and  ion-atom  rearrangement  reactions  may  occur.  Tlie  reaction  rate 

coefficients  for  most  of  these  reactions  are  not  know.  IVe  selected  the 

following  two  near  resonance  charge  exchange  processes  from  a large  set 
1'’ 

given  elsewhere 

+ 0(^D)  - .N(^D)  + 0'"(^S)  (3) 

+ 0(^P)  - + 0‘"(‘^S)  (4) 

and  assumed  a reaction  rate  coefficient  of  5 X 10  cm^/sec  for  each  one. 

This  assumed  rate  coefficient  is  small  compared  with  near  resonance  charge 
exchange  processes  which  are  of  the  order  of  ~ 3 X 10  cm'^/sec.  To  assess 
the  importance  of  reactions  (3)  and  (4),  a given  initial  condition,  Case  II, 
was  used  in  sensitivity  calculations  and  the  results  of  certain  species 
concentrations  are  shown  in  Figures  15-15.  In  Figure  13,  the  electron  density 
calculations  with  and  without  reactions  (3)  and  (4)  are  shown  where  the  lower 
curve  presents  the  electron  density  with  the  inclusion  of  reactions  (3)  and 
(4).  It  is  apparent  that  the  divergence  in  the  electron  density  calculated 
with  and  without  reactions  (3)  and  (4),  increases  with  time.  The  fast  decay 
of  electrons  (lower  curves  of  Figs.  3 and  4)  is  due  to  the  conversion  of 
into  0*,  which  forms  .\0^  in  its  reaction  with  N^.  The  corresponding  SO, 

SO*  and  :N'(^D)  densities  are  shown,  with  and  without  reactions  (3)  and  (4), 
in  Figs.  14  and  15.  Clearly,  if  one  raises  the  rate  coefficient  to  higher 
and  perhaps  more  reasonable  values,  the  divergence  in  densities  would  be 
much  higher.  However,  even  with  a rate  coefficient  as  small  as  5 X 10 
cm'^/sec,  one  can  assess  that  reactions  (3)  and  (4)  are  important  at  late- 
times  for  communications  and  infrared  emission. 
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Fig.  7 - The  electron  density  for  Case  I as  a function  of 
time  using  the  NO+  dissociative  recombination  rate  coef- 
ficients of  Fig.  6.  The  notations  0,  W and  H stand  for 
0 (old  Ref.  5),  W (Walls  auid  Dunn  Ref.  10)  and  H (Huang, 
et  al.  Ref.  9),  respectively.  The  temperature  dependenc 
of  these  rates  are  0 ~ Tg"^*  2,  W ~ Tg"®*  pj  ~ Tg"® 

respectively. 


~ T ”0. 83 


and  H 


Fig.  8 - The  NO"*"  density  for  Case  I using  rate  coefficients 
of  Fig.  6.  The  notations  0,  W and  H are  described  in  the 
caption  of  Fig.  7. 
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Fig,  11  - The  electron  density  for  Case  V using  the 
dissociative  recombination  rates  of  Fig.  6 
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TIME  (SEC) 

Fig.  13  - The  electron  density  as  a function  of  time  for 
Case  n with  and  without  reactions  3 and  4 (see  text).  1 
lower  curve  is  with  reactions  3-  and  4 being  included  in 
the  calculation. 
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TIME  (SEC) 

Fig.  14  - The  NO  and  NO"^  densities  for  Case  II  with  and 
without  reactions  3 and  4 (see  text).  The  upper  curves 
are  with  reactions  3 and  4 being  included  in  the  calculation 
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54^IES 

Offutt  AFB,  NB  68115 

ATTN:  RDPO  LT.  Alan  B.  Merrill 

Los  Alamos  Scientific  Laboratory 

P . 0 . Box  1665 

Los  Alamos,  NM  875^- 

ATTN:  DOC  COM  for  R.  F.  Taschek 
ATTN:  DOC  CON  for  Milton  Peek 
ATTN:  CCC  CON  for  Eric  LindmaJi 
Sandia  Laboratories 
P.  0.  Box  5800 
Albuquerque,  NM  87115 

ATTN:  DOC  CON  for  A.  Dean  Thronb rough 

ATITT;  DOC  CON  for  W.  D.  Brown 

ATTN:  DOC  CON  for  D.  A.  Dahlgren,  ORG  1722 

ATTN:  DOC  CON  for  J.  ?.  Martin,  CRG  1732 

University  of  California 
Lawrence  Livermore  Laboratory 
P.0,  Box  8c8 
Livermore,  CA  9^550 

ATTN:  Tech  Info  Dept  L-5 

Department  of  Commerce 
National  Oceanic  & Atmospheric  Admin, 
Environmental  Research  Laboratories 
Boulder,  CO  80302 

ATTN:  Joseph  H,  Pope 
ATTN:  C,  L,  Rufenach 

Department  of  Commerce 
Office  for  Telecommunications 
Institute  for  Telecom  Science 
Boulder,  CO  60302 

ATTN:  Glenn  Falcon 
ATTIh:  G,  Reed 

AXIN.;  L,  A.  Berr>- 
AiiN:  villiam  F.  Utlaut 

Department  of  Transportation 
Transportation  Rsch . System  Center 
Kendall  Square 
Cambridge,  MA  02lh2 

ATTN:  TER  G.  Harowles 

N.ASA 

Goddard  Space  Flight  Center 
Greenbelt,  Md  20771 

ATTN:  CODE  750  T,  Golden 
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NASA 

600  Independence  Ave . , S . W . 
Washington,  D.  C.  203^6 
ATTN:  M.  Dub in 

Aerodyne  Research,  Inc. 

Tech/Ops  Building 
20  South  Avenue 
Burlington,  MA  OI805 
ATTN:  M.  Camac 
ATTN:  F . Bien 

Aerospace  Corporation 
P.  0.  Box  92957 
Los  Angeles,  CA  90009 
ATTN:  T.  M.  Salmi 
ATTN:  S.  P.  Bower 
ATTN:  V.  Josephson 
ATTN:  SMFA  for  PW^W 
ATTN:  R.  Grove 
ATTN:  R.  D.  Rawcliffe 
ATTN:  T.  Taylor 
ATTN:  Harris  Mayer 
ATTN:  D.  C.  Cartwright 

Analytical  Systems  Corporation 
25  Ray  Avenue 
Burlington,  MA  OI803 

ATTN:  Radio  Sciences 

Avco-Everett  Research  Laboratory,  Inc. 
2385  Revere  Beach  Parkway 
Everett,  MA  021^9 

ATTN:  Richard  M.  Patrick 

Boeing  Company,  The 
P . 0 . Box  3707 
Seattle,  WA  9812^4- 

ATTN:  D.  Murray 
ATTN:  Glen  Keister 

Brown  Engineering  Company,  Inc. 
Cummings  Research  Park 
Huntsville,  AL  35807 

ATTN:  David  Lambert  MS  I8 

California  at  San  Diego,  Univ.  of 
Building  5OO  Mather  Campus 
3172  Miramar  Road 
La  Jolla,  CA  92037 

ATTN:  Henry  G.  Booker 
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Calspan 
P . 0.  Box  235 
Buffalo,  N.  Y.  14221 

ATT1»:  Romeo  A.  Deliberis 

Computer  Sciences  Corporation 

P.  0.  Box  550 

6565  Arlington  Blvd. 

Falls  Church,  VA  22046 
ATTN:  H.  Blank 
ATTN:  Barbara  F.  Adams 

Comsat  Laboratories 
P.  0.  Box  115 
Clarksburg,  Md.  20754 
ATTN:  R.  R-  Taur 

Cornell  University 

Department  of  Electrical  Engineering 
Ithaca,  N.  Y.  14850 

ATTN:  D.  T.  Farley,  Jr. 

ESL,  Inc. 

495  Java  Drive 
Sunnyvale,  CA  ^^102 

ATTN:  J.  Roberts 
ATTN:  V.  L.  Mower 
ATTN:  James  Marshall 
ATTN:  R.  K.  Stevens 

General  Electric  Company 
Tempo-Center  for  Advanced  Studies 
816  State  Street 
Santa  Barbara,  CA  93102 
ATTN:  Don  Chandler 
ATTi;:  DASIAC 
ATTN;  Tim  Stephens 


General  Electric  Company 
P.  0.  Box  1122 
Syracuse,  N.  Y.  15201 

ATTN;  F.  A.  Reibert 

General  Research  Corporation 
P.  0.  Box  5537 
Santa  Barbara,  CA  95^05 
ATTN:  John  Ise,  Jr. 
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Geophysical  Institute 
University  of  Alaska 
Fairbanks,  AK  99701 

ATTN : Technical  Library 
ATTN:  Neil  Brown 
ATTN:  T.  N.  Davis 

GTE  Sylvania,  Inc. 

189  B Street 

Needham  Heights,  MA  02191- 
ATTN:  Marshall  Cross 

HRB -SINGER,  Inc. 

Science  Park,  Science  Park  Road 

P.  0.  Box  60 

State  College,  PA  I68OI 

ATTN : Larry  Feathers 

Honeywell  Incorporated 
Radiation  Center 
2 Forbes  Road 
Lexington,  MA  02175 

ATTN:  W.  Williamson 

Illinois,  University  of 
Department  of  Electrical  Engineering 
Urbana,  IL  6I8OI 

ATTN:  K.  C.  Yeh 

Institute  for  Defense  Analyses 
loo  Army-Navy  Drive 
Arlington,  VA  22202 

ATTN : Ernest  Bauer 
ATTN:  Hans  Wolfhard 
ATTN;  J.  M.  Aein 
ATTN:  Joel  Bengston 

Inti  Tel  & Telegraph  Corporation 
500  Washington  Avenue 
Nutley,  N.  J.  O7IIO 

ATTN:  Technical  Library 

ITT  Electro-Physics  Laboratories,  Inc. 
91IO  Old  Annapolis  Road 
Columbus,  Md.  2101-3 

ATTN;  John  M.  Kelso 
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Johns  Hopkins  UniversiCy 
Applied  Physics  Laboratory 
c621  Georgia  Avenue 
Silver  Spring,  MD  20910 

ATTN:  Docunent  Librarian 

Lockheed  Missiles  & Space  Co.,  Inc. 

P.  0.  Box 

Sunnyvale,  CA  9^23 

ATTN:  Dept.  60-12 

Lockheed  Missiles  and  Space  Company 
5251  Hanover  Street 
Palo  Alto,  CA  9^30h- 

ATTN:  Billy  M.  McComac,  Dept  52-1^ 
ATTN;  Martin  Walt,  Dept  52-10 
ATTN:  Richard  G.  Johnson,  Dept  52-12 
ATTN;  JOHN  CLADIS 


MIT  Lincoln  Laboratory 
P.  0.  Box  73 
Lexington,  MA  02173 

ATTN:  Mr.  Walden,  XII3 

ATTN:  D.  Clark 

ATTN:  James  H.  Pannell,  L-2h6 

ATTN;  Lib  A-0S2  for  David  M.  Towle 

Martin  Marietta  Corporation 
Denver  Distribution 
P.  0.  Box  179 
Denver,  CO  60201 

ATTN:  Special  Projects  Program  248 

Max-.^ell  Laboratories,  Inc. 

92I4  Balboa  Avenue 
San  Diego,  CA  92125 

ATTl'J:  A.  J,  Shannon 
ATTN:  V.  Fargo 
ATTN;  A.  N.  Rostocker 

McDonnell  Douglas  Corporation 
5301  Bolsa  Avenue 
Huntington  Beach,  CA  92637 
ATTi;:  J.  Moule 
Am:;  N.  Harris 
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Mission  Research  Corporation 
735  Stats  Street 
Santa  Barbara,  CA  93101 
ATTI'I:  R.  Hendrick 
ATTH:  Conrad  L.  Longaire 
ATTIT:  Ralph  Kilb 
A.TTTI;  R.  E.  Rosenthal 
ATTIT:  R.  Bogusch 
ATTiT : David  Sc  wle 
ATTH:  M.  Scheibe 
ATTH:  P.  Fischer 
Mitre  Corporation,  The 
Route  62  and  Middlesex  Turnpike 
?.  0.  Box  208 
Bedford,  liA  01730 

ATTN:  Chief  Scientist  W.  Sen 
ATTN:  S.  A.  Morin  h/S 
ATTN:  C.  Ilirding 


North  Carolina  State  Univ  At  Raleigh 
Raleigh,  N.  C.  27507 

ATTN:  SEC  Officer  for  Walter  A. 

Pacific-Sierra  Research  Corp. 

IU56  Cloverfield  Blvd. 

Santa  Monica,  CA  9040U 

ATTN:  E.  C.  Field,  Jr. 

Philco-Ford  Corporation 
Western  Development  Laboratories  Div 
3939  Fabian  Way 
Palo  Alto,  CA  9I303 

ATTN:  J.  T.  Mattingley  MS  X22 

Photometries,  Inc. 

IU2  ITarrett  Road 
Lexington,  I'lA  02173 

ATTN:  Irving  J.  Kofsky 

Mitre  Corporation,  The 
Westgate  Research  Park 
1^20  Dolley  Madison  Blvd. 

McLean,  VA  22101 

A.TTN:  Allen  Schneider 


Flood 
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Physical  Dynamics j Inc. 

?.  0.  Box  1059 
Berkeley,  CA  9^701 

ATIOJ:  Joseph  B.  Worlcnan 

Physical  Sciences,  Inc. 

607  North  Avenue,  Door  iS 
Wakefield,  MA  0l880 
ATTN:  ’vray 

R & D Associates 

P.  0.  Box  55cO 

Santa  Monica,  CA  90^3 

ATTN;  Robert  E.  Lelevier 
ATTN:  Forest  Gilmore 
ATTN:  Richard  Latter 
ATTN:  William  B.  Wright,  Jr. 

Pv  & D Associates 
iSlp  N.  Ft.  Myer  Drive 
11th  Floor 
Arlington,  VA  2209 

ATTN:  Herbert  J.  Mitchell 

Rand  Corporation,  The 
1700  Main  Street 
Santa  Monica,  CA  $0406 
ATTN:  Cullen  Crain 


Science  Applications,  Inc. 

P.  0.  Box  2551 
La  Jolla,  CA  92058 

A-TTN:  Daniel  A.  Hamlin 
ATTN:  D.  Sachs 
ATT:;:  E.  A.  Straker 

Space  Data  Corporation 
1531  South  26th  Street 
Phoenix,  AZ  8503^ 

ATTN:  Edward  F.  Allen 
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Star.ford  Research  Institute 
333  Ravenswood  Avenue 
Park,  CA  9^025 
ATH'T:  M.  Earon 
ATTII:  L.  L.  Cobb 
A?x?I:  Walter  G.  Chestnut 
ATT?!:  David  A.  Johnson 
ATTIT:  Charles  L.  Rino 
ATT?!:'  E.  J.  Fremouw 
ATT?!:  Ray  L.  Leadabrand 
ATTN:  Donald  Neilson 

Stanford  Research  institute 
3C6  Wynn  Drive,  N.  W. 

Huntsville,  AL  35805 

ATT?!:  Dale  H.  Davis 

Technology  International  Corporation 
75  Wiggins  Avenue 
Bedford,  ?!A  01730 

ATTN:  W.  P.  Boquist 

TR'W  Systems  Group 
Cne  Space  Park 
Redondo  Beach,  CA  90278 
ATT?!:  P.  H.  Katsos 
ATT?!:  J.  W.  Lowery 

Utah  State  University 
Logan,  UT  8U32I 

ATT?!:  C.  .Wyatt 
ATTN:  D.  Burt 
ATTIT:  Kay  3al;er 
ATTN : Doran  BaJcer 

7isid3me,  Inc. 

19  Third  Avenue 

North  West  Industrial  Park 

Burlington,  MA  OI803 

ATT?!:  William  Reidy 
ATTN:  Oscar  Manley 
ATT?!:  J.  W.  Carpenter 
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